Cosmic spherical void via coarse-graining and averaging non-spherical structures. 
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Abstract 



" |~ Inhomogeneous cosmological models are able to fit cosmological observations without dark energy under the assumption that we 
live close to the "center" of a very large-scale under-dense region. Most studies fitting observations by means of inhomogeneities 
also assume spherical symmetry, and thus being at (or very near) the center may imply being located at a very special and unlikely 
observation point. We argue that such spherical voids should be treated only as a gross first approximation to configurations that 
C~\ follow from a suitable smoothing out of the non-spherical part of the inhomogeneities on angular scales. In this Letter we present 
D a toy construction that supports the above statement. The construction uses parts of the Szekeres model, which is inhomogeneous 
l~H and anisotropic thus it also addresses the limitations of spherical inhomogeneities. By using the thin-shell approximation (which 
means that the Israel-Darmois continuity conditions are not fulfilled between the shells) we construct a model of evolving cos- 
mic structures, containing several elongated supercluster-like structures with underdense regions between them, which altogether 
provides a reasonable coarse-grained description of cosmic structures. While this configuration is not spherically symmetric, its 
proper volume average yields a spherical void profile of 250 Mpc that roughly agrees with observations. Also, by considering a 
^•■^ non-spherical inhomogeneity, the definition of a "center" location becomes more nuanced, and thus the constraints placed by fitting 
observations on our position with respect to this location become less restrictive. 
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Cosmological models allowing for non-trivial inhomo- 
geneities have become a popular tool to analyze cosmological 
observations without the need of introducing an elusive dark 
energy source (for a review on a subject and explicit exam- 
ples the reader is referred to Ref. [Q]]). Within this approach, 
the preferred configurations are Gpc-scale cosmic void mod- 
els based on the spherically symmetric Lemaitre-Tolman (LT) 
models , under the assumption that we live close to a cen- 
ter of a cosmic depression of radius around 1-3 Gpc J3, [5I [(J . 
These configurations are often criticized on the grounds that 
they violate the Copernican principle, since compliance with 
the cosmic microwave background (CMB) constraints implies 
that only one such Gpc structure is allowed and the observer 
cannot be further away from the origin than ~ 50 Mpc [7] Q. 
However, as suggested recently by Alexander et al. [8] (see 
also |@1), a void of radius 250 Mpc is sufficient to explain the 
supernova observations, the power spectrum of the cosmic mi- 
crowave background and is also consistent with Big Bang Nu- 
cleosynthesis, or Baryon Acoustic Oscillations. For void struc- 
tures of this size the Copernican Principle is not violated - there 
can be many such structures as the upper size to violate the 



CMB constrains is 300 Mpc |T(1 U]. Also, restricting our po- 
sition to be within 50 Mpc from the origin of a 250 Mpc void is 
a less stringent constraint. However, one may argue that voids 
of radius 250 Mpc are not observed in the filamentary structure 
characterizing our Local Universe. 

In this Letter we address this issue by showing that these 
rather artificial spherical void structures need not exist in its 
pure form. Instead, they approximate configurations that can 
emerge after coarse-graining and averaging a sufficiently large 
scale region of a realistic lumpy Universe in which the density 
distribution is far from spherical. We consider for this purpose a 
non-spherical inhomogeneous and anisotropic Szekeres model, 
prescribing its free functions by means of a thin-shell approxi- 
mation [ 1( 



111 LL2L LL2I], leading to a reasonable coarse-grained 
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'The figure obtained in Ref. 01 is 15 Mpc, but as shown in |6] fluctuations 
of Hubble rate needed to fit observations are Ah x 0.08 which means that for 
d < 55 Mpc the dipole is less than 3 mK. 



description of realistic cosmic structures. Initial conditions are 
defined at the last scattering surfaces to show that the model 
can evolve from small early universe initial fluctuations, hence 
we achieve consistency with models of structure formation. We 
show that averaging this inhomogeneous non-spherical config- 
uration leads to a cosmic void that is qualitatively similar to the 
spherical models discussed by Alexander et al. (8Q. 

It is important to emphasize that the lack of spherical symme- 
try in the Szekeres model removes the unique invariant nature 
of the center location of models with this symmetry. This is 
a relevant feature of the model, as our being sufficiently near 
this center is a strong constraint that the fitting of observations 
place on spherical LT models. As we argue further below, this 
constraint may be less restrictive in a non-spherical Szekeres 
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model. 



2. Szekeres model 



2.1. Einstein's equations 

The metric of the Szekeres model takes the following form 

Hi 



ds z = At 



2 (g>'-t&g7g) 2 dj2 



o 2 



-(dx 2 + dy 2 ), 



(1) 



where <S> = <5(f, r), the prime denotes the partial derivative with 
respect to r, O' = d$>/dr, and 



£= - 



+ e 



(2) 



while k(r), S (r), P(r), Q(r) are arbitrary functions; e is a con- 
stant: the e = — 1 case is often called the quasi-hyperbolic Szek- 
eres model, e = quasi-plane, and e = 1 quasi-s phe rical (for a 



detailed discussion on these models see II15L llq, 11711 ') - in this 



Letter we only consider the quasispherical case. The coordinate 
system in which the metric has the form ([T} can be interpreted 
as a stereographic projection of polar coordinates 1171 12511 . For 
the quasispherical case the transformation is of the following 
form: 

{x-P,y- Q) = {scot(|)cos(0),Scot(|)sin(, 



Then, using the (8, <f>) coordinates we can rewrite G'/G as [17 
25] 



G' S, r cos 8 + sin 9 (P, r cos <p + Q, r sin <p) 
~G = S ■ 



(3) 



However, under the above transformation the metric becomes 
non-diagonal |2(| 21]. Thus, for some applications it is more 
convenient to use the (f , r, x, y) coordinates - see Sec. [3] Notice 
that we use the (f, r, 8, <p) coordinates in Sec. 12.41 

Einstein's equations for a dust source associated with ([T])-© 
reduce to 



O 2 = -k(r) + 



2M(r) 
<t> ' 



(4) 



where M(r), is an arbitrary function related to the density p via: 
2 AT - 6MG'/& 



k P = 



<t> 2 (<6' -<&£'/£)' 



(5) 



where k = 8ttG. Whenever O' = <t>£'/£ and M' + 3MG'/G 
hold a shell crossing singularity occurs, which (in a quasi- 
spherical model) may occur along a circle, or, in exceptional 
cases, at a single point, and not at a whole surface in the (t, r) 
plane, as is the case in LT models 11171 12511 . 

As in a spherical LT model, the bang time function follows 
from © 



d<D 



= t - t B (r). 



(6) 



k + 2M/<D 



2.2. Set-up 

While the previous equations indicate that the Szekeres 
model is specified by 6 functions, choosing appropriate coor- 
dinates eliminates one of the independent functions. Thus, we 
must provide 5 functions as initial conditions to specify the 
model. In particular, we will specify the functions S,P, Q,ts, 
and M. 

In order to achieve with a Szekeres model the most realistic 
possible description of cosmic structures and structure forma- 
tion, we specify below suitable forms for these 5 functions and 
define our model at the last scattering instant, so that it evolves 
from small initial fluctuations up to the present cosmic time. 
The algorithm that we use in the calculations can be defined as 
follows: 

1 . The radial coordinate is chosen to be the areal radius at the 
last scattering instant f - <£(r, f,). However, to simplify 
the notation we will omit the bar and denote the new radial 
coordinate by r. 

2. The chosen asymptotic cosmic background is an open 
Friedman mode^ i.e. Q„, = 0.3 and A = 0. The back- 
ground density is then given by 



Pb 



= O, 



3H~ 

x Pc ,. = 0.3x— ^(1+z) 3 , 
isnG 



(7) 



where the Hubble constant is Hq = 70 km s 1 Mpc 1 . 

3. The age of the universe is given by (|6]). Notice that inho- 
mogeneities may affect the age relation via the bang time 
function tg. However, since we define our model at the last 
scattering when the Universe is expected to be very close 
to homogeneity, we assume that fg = 0. As a consequence, 
the age of the Universe is everywhere the same (as in the 
homogeneous background Friedmann model) and is equal 
to t t = 471, 509.5 yean0 (see d for details). 

4. The function M{r) is given by 



M(r) = 4n- 



f 
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p b {\+5p)r dr, 



where 6p = -0.005e-W 100 > 2 + 0.0008e-K f - 50 )/ 35 ] 2 + 
0.0005e-K f - I15) / 6 °] 2 + 0.0002e-K f - 14( W 55 l\ and I = r/ 1 
kpc. 

5. The function k(r) can be calculated from ([6j. 

6. The functions Q, P, and S are defined as follows 
s = 1 =>S' = 0, 

D= 1.05(1 + r)- a »e-°- 004r , 



2 The standard conventional analysis of cosmological observations appears 
to imply a spatially flat background. However, when the assumption of homo- 
geneity is relaxed spatial flatness is no longer required. For specific examples 
and a comprehensive discussion on this issue the reader is kindly referred to 

da. 

3 The age of the universe at the last scattering instant is below 400, 000 yr. 
This is due to the presence of radiation that is not negligible at the decoupling 
instant. However, the Szekeres model only describes the pressureless config- 
uration (dust). That is why when calculating the initial instant we neglect the 
pressure and only consider dust. If radiation is neglected then the age of the 
universe at the decoupling instant is larger than 400, 000 yr. 
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Q =£>, P' = 0for£s^27, 

Q' = -2), P' = for 27 < f < 35, 

Q' = 0, P' = -2) for 35 < € < 41, 

Q' =0, P' = 2) for 41 < ^ ^ 51.5, 

Q! = 0.882), P' = -0.52) for 51.5 <€^6l, 

Q' = 0.712), P' = 0.712) for 61 < I < 69, 

Q' = 0, P' = -2) for 69 < £ < 77, 

G' = -2), P' = for 77 < £ < 86.5, 

Q = 0.1 AD, P' = -0.7 '42) for 86.5 < £ sC 96, 

Q = D, P' = 2) for 96 < £ < 102, 

Q' = -2), P' = for 102 < i < 1 15, 

Q' = 2), P' = for 115 < t < 129, 

g' =0, P' = -2)for£> 129. 

7. Once the model is specified, its evolution is calculated 
from eq. ©. 

8. Density distribution at the current instant is then evaluated 
from ©. 
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As can be seen from the expressions above, some of the 
functions are not continuous, hence the Israel-Darmois match- 
ing conditions for spacetimes are not fulfilled. Therefore, this 
model should be regarded as a thin-shell app roximation. This 
approximation, in its original form II 1(1 ll U Il2l 1 1 311 - consists 
of: (i) two Friedmann regions and (ii) transitions between them 
modeled by thin shells with negligible effects. Here we apply a 
similar approach, but instead of Friedmann regions we consider 
Szekeres regions. This approach is justified by the fact that in 
the Szekeres models each surface of constant r (a 2-sphere in 
our case) evolves as an independent Friedmann model - as seen 
from eq. <j4j each shell can be considered as a Friedmann model 
whose evolution is given by its Friedmann equation ©, a dif- 
ferent Friedmann equation for each sphere. The only restriction 
is that the model must be free from shell crossings - an inner 
sphere cannot expand faster than an outer sphere. This has been 
taken into account and the above defined model evolves without 
shell crossings from the initial instant (last scattering) until the 
present moment. 

The density distribution for our model is depicted in Fig. Q] 
(we remark that this model does not admit Killing vectors). We 
are using intuitive Cartesian coordinates (see Ref. 

MM for 

the corresponding transformation and description of these co- 
ordinates). It should now be clear why we have selected the 
above set of functions to define the Szekeres model. As seen 
in Fig. [1] we have all together 13 different sets with the func- 
tions {M, tb, Q,P,S} and so we have 13 different overdensities. 
By changing the form of the above functions we can change the 
pos ition, size, and the amplitude of the overdensities (see Ref. 
B20II21I1 for a detailed discussion). It is important to remark that 
the very special choice of free functions is not an act of "fine 
tuning" (in the sense of a tricky manipulation of parameters), 
but an effort to approximate the known Cosmography as best 
as possible. Hence, if new data would modify or improve the 
known density distribution it would be straightforward to mod- 
ify these functions to approximate this data, which would then 
lead to a structure that is different from that presented in Fig. Q] 



Figure 1: The present-day color-coded density distribution p/po (where pg 
is density of the homogeneous background model). Brighter colors indicate 
a high-density region, darker low-density region. The presented surface is a 
horizontal cross-section, i.e. 8 = n/2. For a discussion on how realistic this 
structure is see Sec. 12.31 



2.3. How realistic this model is ? 

As shown in Fig. [1] the model under consideration contains 
structures such as voids and elongated supercluster-like over- 
densities. It has large overdensities around ~ 200 Mpc ( on the 
left in Fig. []]) that compensate the underdense regions and al- 
low the model to be practically homogeneous at r > 300 Mpc. 
In the real universe we observe very massive matter concentra- 
tions - the Shapley Concentration roughly at the distance of 200 
Mpc, or the Great Sloan Wall at the distance of 250-300 Mpc. 
In the opposite direction on the sky we find the Pisces-Cetus 
and Horologium-Reticulum, which are massive matter concen- 
trations located at a similar distance. We refer the reader to Fig. 
44 of Ref. B23I1 . which provides a density map of the Local 
Universe reconstructed from the 2dF Galaxy Redshift Survey 
Survey using Delaunay Tessellation Field Estimatofl Also, the 
inner void seen in Fig. Q]is consistent with what is observed in 
the Local Universe - it appears that our Local Group is not lo- 
cated in a very dense region of the Universe, rather it is located 
in a less dense region surrounded by large overdensities like the 
Great Attractor on one side and the Perseus-Piscis supercluster 
on the other side. Both are located at around 50 Mpc — see 
Fig. 19 of l24ll that provides the density reconstruction of the 
Local Universe using the POTENT analysis. 

While still far from a perfect "realistic" description, the 
model displayed in Fig. [T]exhibits the main features of our local 
Universe. It should be therefore treated as a "coarse-grained" 
approximation to study local cosmic dynamics by means of a 



4 This figure is also available at http://en.wikipedia.Org/wiki/File:2dfdtfe.gif 
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suitable exact solution of Einstein's equations. Such approxi- 
mation is, evidently, far less idealized than the gross one that 
follows from spherically symmetric LT models. 

2.4. Position of the "center" 

Since the model we are considering is not spherically sym- 
metric, there is no invariant and unique characterization of a 
center worldline. Instead, for every 2-sphere corresponding to 
a fixed value of r, we have (at least) two locations that can be 
considered appropriate generalizations of the spherically sym- 
metric center: the worldline marked by r — where the shear 
tensor vanishes, which defines a locally isotropic observer (cf. 
eq (16.29) of Ref. [25]), and the "geometric" center of the 2- 
sphere whose surface area is An<& 2 . 

Since the 2-spheres of constant r in a quasi-spherical Szek- 
eres model are non-concentric, their geometric center does not 
coincide with the the point of vanishing shear r = - this is 
schematically presented in Fig. [2] As a result, the distance from 
the origin to the surface of the sphere depends on the direction 
marked by the angles [8, (/>): 



Sir, 6, 
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d>' -$£'/£ 
dr — , 



(8) 



and thus, the displacement A between the origin and the geo- 
metric center of a sphere marked by comoving radius r is given 
by 




Figure 2: Schematic representation of the quasispherical Szekeres model - 
different surfaces of constant t and r (spheres) have different centers, which 
do not coincide with the origin marked by r = (denoted by a black dot). 
The center of the larger sphere (geometric center) is depicted by a cross. The 
maximal distance from the origin to the surface of the sphere is 5 max and the 
minimal is S mm . The distance between the geometric center and the origin is 
denoted by A. 



A = 



where (5 max = max (5), 6^ m = min(5). As seen from (0, the 
maximal and minimal value of £' /£ for our model (where S ' = 
0) corresponds to 6 — n/2. The distance, 6, as a function of 
for 4 cases [(1) present-day areal radius is Q)(to,r) = 100 
Mpc, (2) <D(f , r) = 150 Mpc, (3) O(f , r) = 200 Mpc, and (4) 
<l>(fo, r) = 250 Mpc] is depicted by Fig. [3] As seen, for example, 
for a sphere whose present-day area radius is <t> = 100 Mpc the 
model under consideration yields a displacement of A = 36 
Mpc towards <p w 80° direction. While for O = 250 Mpc we 
have A = 62 Mpc towards ip w 120°. 

Different 2-spheres of constant r have their centers displaced 
(with respect to the coordinate origin where the shear vanishes) 
with different values and towards different directions. Hence, 
there is no uniquely defined geometric center - a center of one 
sphere does not coincide with a center of another sphere and 
moreover it also does not coincide with the point where shear 
vanishes. Since the center in spherical LT models is also a lo- 
cally isotropic observer (where shear vanishes), and fitting ob- 
servations in these models restricts our cosmic location to be 
within a given maximal separation from this observer, then it 
is reasonable to expect that similar restrictions should emerge 
in a Szekeres model given in terms of maximal separation from 
the local isotropic observer at the coordinate origin where shear 
vanishes. Therefore, the displacement of the geometric center 
from this origin would make our location in a Szekeres model 
less special and improbable than in spherically symmetric mod- 
els where both center locations (local isotropic observer and ge- 
ometric center) necessarily coincide. 




Figure 3: Proper distance, S, given by (8) as a function of (/> for 4 different cases: 
tD(f , r) = 250 Mpc (the dash-dotted line), <J>(f , r) = 200 Mpc (the dotted line), 
<D(fo, r) = 150 Mpc (the dashed line), <D(fo, r) = 100 Mpc (the solid line). For 
all cases 6 = n/2. 
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3. Averaging 

As shown in Ref. I22I1 . the proper 3-dimensional volume 
in space slices orthogonal to the 4-velocity (t = constant) in a 
Szekeres model is 



V© = Jdrjdxj dyV=S 

—co —00 

= 4n dr = 

J vT^I 



(9) 



and thus, the proper volume averaged density is spherically 
symmetric (i.e. independent of x and y), even if the density 
itself is far from a spherical distribution: 



DO DO 



() —CO —CO 



<PX>©) = ^- J dr j dx j ' dy^Pgp 

1 r 2AT 
I dr 

kRd J Vl -k 



(10) 



This spherical volume-averaged density distribution, (p)(rx)), 
evaluated as a function of ro, is displayed by Fig. [4] while its 
profile along the radial direction is presented in Fig. [5] As can 
be seen from these figures, the "angular" (i.e. x,y) dependence 
of an inhomogeneity (which is a highly non-spherical coarse 
grained density distribution) has been smoothed out by the av- 
eraging process, resulting into an averaged distribution that is 
equivalent to a spherical cosmic void whose radius is approxi- 
mately 250 Mpc (as in Ref. H). 

It is important to remark (and clarify) several issues. First, 
the field source is the local non-spherical density, p, not its 
spherically symmetric average, (p). Hence, the former (and 
not the latter) must be considered in the study of the general 
relativistic dynamics of the inhomogeneity, as well as in all 
computations related to observations (based on null geodesies). 
Second, considering the dynamics of (p) by means of the sug- 
gested averaging procedures for General Relativity ll26ll is be- 
yond the scope of this Letter. Instead, since (p) follows from 
a proper volume averaging procedure that is consistent for co- 
variant scalars, we simply regard this averaged density as an 
approximation that conveys useful non-local information on p, 
even if it cannot account for the latter's description of the intri- 
cacies of local detail. 

While the proper volume average (p) in quasi-spherical 
Szekeres is necessarily a spherical distribution [22], it is not 
obvious (and not possible to know beforehand) the form of its 
radial density profile, though some sort of void profile should be 
expected given the observed void dominance in cosmic struc- 
ture. In this context, we regard as an important result the 
fact that the simple shallow void profile displayed in figure 
[5] emerged by averaging the specific non-spherical p that we 
proposed, which, as mentioned before, provides a reasonable 
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Figure 4: Volume average density distribution (p)(rj})/po (where po is the back- 
ground density). Brighter colors indicate a high-density region, darker low- 
density region. The displayed surface is a horizontal cross-section, i.e. S = n/2. 
As follows from 11 OK all "angular" dependence of inhomogeneities has been 
smoothed out. 



course grained description of realistic cosmic structure. The 
fact that the void profile of figure [5] (and not any other more 
complicated void-like or clump/void mixed profile) resulted is 
an encouraging signal, as it provides a theoretical connection 
with spherical void models (as those studied in 1 8] whose pro- 
files resemble that of figure 0. We also have a concrete ex- 
ample suggesting that spherical models, as idealized approxi- 
mations, may be analogous to the average (first integral mo- 
mentum approximation) of well motivated (and less idealized) 
non-spherical models. As a consequence, the use of a Szekeres 
model seems to suggest that results obtained by means of spher- 
ical LT models may be robust: while local non-spherical infor- 
mation could still provide important refinements, and is needed 
for computations involving null geodesies (specially when fit- 
ting CMB constraints), it is likely that basic bottom line infor- 
mation is already contained in the spherical voids constructed 
with LT models. 



4. Conclusions 

We have presented in this Letter an explicit example of how 
a non-spherical construction using parts of an exact solution 
of Einstein's equations can provide a reasonable coarse-grained 
description of a present day cosmic structure that has evolved 
from small initial fluctuations at the last scattering surface. This 
coarse grinned density distribution produces (after averaging) 
a cosmic void profile that is analogous to density profiles in 
spherical voids of radius ~ 250 Mpc. As mentioned before, 
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Figure 5: Density profile after averaging (p)(/"j))/po (where po is the back- 
ground density). Since (p)(r£>) in ilO\ is independent of x and y, its depen- 
dence on r has the form of a cosmic void whose radius is approximately 250 
Mpc (analogous to the spherical voids of reference 18|]). 



a void of at least this size is required to explain cosmological 
observations without the need for dark energy. 

We emphasize that we do not claim that the coarse-grained 
description given by our model is completely realistic, though 
it contains essential features and thus represents a significant 
improvement over the description of observed cosmic struc- 
ture in spherical LT models used in the literature. Our model 
uses the thin shell approximation which means that the Israel- 
Darmois conditions are not fulfilled between the shells. Also, 
the model becomes practically homogeneous at a distance of 
~ 350 Mpc, and thus it only describes one small region of our 
Universe, which should be composed of many other similarly 
sized structures. Such a construction would be consistent with 
supernova observations and the power spectrum of the cosmic 
microwave background, and would not violate the Copernican 
Principle and the constrains form the Rees-Sciama effect - to 
violate the R-S effect eve ry sing le component structure must 
be larger than ~ 300 Mpc flolfTlfl . 

The model we have presented is one among the first attempts 
in using the Szekeres solution as a theoretical and empiric tool 
to study and interpret cosmological observations [27, 28, 9T. l29tl . 
This opens new possibilities for inhomogeneous cosmologies, 
as this is the most general available cosmological exact inhomo- 
geneous and anisotropic solution of Einstein's equations. The 
model provides a more nuanced and much less restrictive de- 
scription of the need to constrain our location with respect to 
a center location. It is also a concrete example that illustrates 
the possibility that a mildly increasing void profile (required by 
observations) can emerge if local structures are coarse-grained 
and then averaged. Of course, notwithstanding these appeal- 
ing features, the model and its assumptions must be subjected 
to hard testing by data from the galaxy redshift surveys, and 
evidently the more comprehensive this data can be the better 
it can be used for this purpose. Unfortunately current surveys 
like 2dF of SDSS do not cover the whole sky and only focus 
on small angular regions of it. However in the near future this 



limitation may be overcome - for example, Sky Mapped aims 
to cover the whole southern sky which will provide sufficient 
data to test possibilities suggested and elaborated in this Letter. 
A more comprehensive and detailed article on the model pro- 
posed here is currently under elaboration and will be submitted 
soon for publication. 
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